Discrete breathers in Fermi-Pasta-Ulam lattices 
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We study the properties of spatially localized and time-periodic excitations - discrete breathers 
- in Fermi-Pasta-Ulam (FPU) chains. We provide a detailed analysis of their spatial profiles and 
stability properties. Especially we demonstrate that the Page mode is linearly stable for symmetric 
FPU potentials. A resonant interaction between a localized and delocalized perturbations causes 
weak but finite strength instabilities for asymmetric FPU potentials. This interaction induces Fano 
resonances for plane waves scattered by the breather. Finally we analyse the interplay between 
energy thresholds for breathers in the presence of strongly asymmetric FPU potentials and the 
corresponding profiles of the low-frequency limit of breather families. 
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The Fermi-Pasta-Ulam (FPU) paradox was ob- 
served fifty years ago. The surprising finding was 
a localization of energy in the reciprocal (/-space 
of a model with discrete translational invariance, 
despite the presence of interaction between ex- 
tended normal modes. Thirty three years later 
Sievers, Takeno and Kisoda reported on the ob- 
servation of energy localization in real space for 
the same class of FPU models, which is as sur- 
prising since these excitations, called discrete 
breathers or intrinsic localized modes, violate the 
underlying discrete translational symmetry of the 
model. The past decade has witnessed a tremen- 
dous progress in the theory and applications of 
discrete breathers, which goes much beyond the 
scope of the original FPU frame. We use the mod- 
ern theory of discrete breathers to investigate the 
properties of these solutions in FPU models, pay- 
ing special attention to the issues of stability, res- 
onances, wave scattering and energy thresholds. 



INTRODUCTION 

The celebrated Fermi-Pasta-Ulam (FPU) model was 
introduced fifty years ago in order to study the process 
of equilibration of energy among normal modes due to 
mode- mode interactions 0. It can be viewed as a toy 
version of a model describing the dynamics of lattice vi- 
brations of perfect crystals. The original FPU model 
reduced space dimension to one, and attributed one de- 
gree of freedom to each lattice site. It is worth to men- 
tion that in the absence of mode-mode interactions, that 
model is often used in solid state physics textbooks to 
explain the basic features of phonons and is then called a 
monoatomic chain Q . Some of these textbooks deal also 
with mode- mode interactions, coining them anharmonic 
corrections (since they appear as anharmonic terms in 
the Taylor expansion of the potential energy of the lat- 



tice with respect to atomic displacements). Notably these 
anharmonic terms are used to explain thermal expansion 
of crystals, among many other features. It is also worth 
to stress that solid state physicists did not seriously ques- 
tion the assumed fact that in thermal equilibrium the en- 
ergy is cquipartitioned among the normal modes - due to 
the bulk of accumulated experimental evidence for that 
fact. The more time- resolved spectroscopical methods 
advance, the more questions arise concerning these top- 
ics when studying lattice dynamics in the presence of 
strong nonlinearities (i.e. mode- mode interactions). Re- 
turning to the FPU model calculations, they revealed 
the surprising fact that at least for some cases (for some 
sets of initial conditions) the evolution of the FPU model 
showed no equipartioning among the normal modes Q]. 
In other words, the FPU paradox seems to consist pri- 
marily of the observation of localization of energy in a 
few normal modes, despite the presence of interaction 
between all modes, which would be capable of distribut- 
ing the energy among all normal modes of the system 

While the modern view on this FPU paradox obser- 
vation is for sure discussed in other contributions to this 
focus issue, we will not dwell on that further. In what fol- 
lows, we will discuss another interesting aspect of the dy- 
namics in the FPU model and its various generalizations. 
This concerns localization in real space. In other words, 
we will show that FPU models allow for solutions of the 
corresponding equations where the energy is not equally 
shared among the local constituents (using the solid state 
physics language, among the different atoms). These so- 
lutions are coined Intrinsic Localized Modes (ILM), or 
Discrete Breathers (DB). Whatever the nomenclature, 
ILMs/DBs turn out to be generic to a much larger class 
of Hamiltonian lattices, and FPU models together with 
their generalizations represent one of the subclasses of 
these lattices (see Fig^for a schematic representation of 
various DBs). ILMs/DBs are time-periodic and spatially 
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where Xi describes the scalar displacement of a particle 
(atom) from its equilibrium position, pi = xi = dxi/dt is 
its conjugated momentum (velocity), I denotes the num- 
ber of the particle, and W{x) is the interaction potential 
between nearest neighbours. The Hamiltonian H is as- 
sumed to take only nonnegative values for small values 
of the displacements and velocities, i.e. the potential W 
and its first derivative W vanish for zero displacements 
W{Q) — W^'(O) — and the second derivative W" is pos- 
itive for small displacements W"{{)) > 0. The Hamilto- 
nian equations of motions xi = dH/dpi , pi = —dH/dxi 
lead to the following set of coupled differential equations: 

ii = -W'{xi-xi^i) + W'{xi+i-xi) . (2) 
Let us expand the function W in the following series 




FIG. 1; A schematic representation of different types of dis- 
crete breathers: (a) acoustic FPU breather; (b) acoustic ro- 
tobreather; (c) optical breather (for details see e.g. Q). 



localized solutions, and exist thanks to the interplay be- 
tween nonlincarity and discreteness |J|. Many studies 
of DBs have been successfully launched, on such topics 
as rigorous existence proofs, dynamical and structural 
stability and computational methods of obtaining DBs 
in classical models as well as their quantum aspects. In 
addition DBs have been detected and studied experimen- 
tally in such different systems as interacting Josephson 
junction systems 0, coupled nonlinear optical waveg- 
uides 6], lattice vibrations in crystals 7], antiferromag- 
netic structures micromechanical cantilever arrays 
91, Bose-Einstein condensates loaded on optical lattices 
lo| . layered high- Tc superconductors DBs are pre- 
dicted also to exist in the dynamics of dusty plasma crys- 
tals It is an equally interesting question of why 
ILMs/DBs have not been properly discussed when study- 
ing the lattice dynamics of crystals until very recently. 
We leave the answer to this question to the experienced 
and educated experts of that field. 



SETTING THE STAGE AND EARLY RESULTS 

We will consider the following class of one-dimensional 
Hamiltonian chains 



pf + W{xi-xi^i) 



(1) 



W{x) 



m=2 



(3) 



If the potential W{x) is symmetric W{x) = W{—x) it 
follows (j)2m+i = for all positive integers m. 

The celebrated FPU models "1| are obtained by choos- 
ing nonzero values for (/)2,3,4 and zeroing all the others 
0m>4 — 0. Thus we may consider the FPU models as low 
amplitude expansions of the more general model class l(T]l. 
FPU models take into account the two first anharmonic 
corrections (^3^4 to the harmonic term (j)2. The case of 
a symmetric potential for FPU models is thus obtained 
by assuming (j)^ = 0. Note that anharmonic terms in the 
Hamiltonian correspond to nonlinear terms in the equa- 
tions of motion here. The equations Q conserve both 
the total energy H ^ as well as the total mechanical 
momentum P = "^iPi- Without loss of generality we 
will consider below the case P = only. 

If we restrict our consideration to very small ampli- 
tudes and velocities only, we may neglect all nonlinear 
terms from the equations of motion, assuming 4>m>2 = 0. 
The solution of the corresponding linear coupled differ- 
ential equations 



xi = 4>2{xi+i + Xl-i - 2xi) 



(4) 



can be written in the form of a superposition of plane 
waves each given by 



Xl{t) — Aq C0s{ujqt — ql + Bq) 



(5) 



where q is the wave number, ujq the plane wave frequency, 
and Aq and Bq are integration constants. The dispersion 
relation (Fig|2) 
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is periodic in q and is characterized by an upper bound 
l^^gl < i^TT = 2\/02. Moreover for small values of q 
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FIG. 2: Dependence of u>q on q for (j)2 = 1- Vertical and 
horizontal lines mark the values q = and ujq = 0. 

\ujq\ ~ Such an acoustic type of dispersion at 

small values of q is intimately connected to the above 
mentioned conservation of the total mechanical momen- 
tum P. Any spatially localized initial excitation on the 
lattice will ultimately disperse due to the g-dependence 
of the group velocity 

v,^ — = ±^,cos-. (7) 

Taking into account weak anharmonic/nonlinear effects 
usually leads to a renormalization of plane wave proper- 
ties such as frequencies and lifetimes. However no princi- 
pal change in the behaviour of a spatially localized initial 
condition is expected from that perspective. 

Some rigorous results exist for the study of stability 
of certain plane waves. Indeed, the corres pon ding mod- 
ulational instability (cf. HI , Q , HI , HI , [H and [13) 
is responsible for the appearance of localized structures. 
However the linearization of the phase space flow around 
the plane wave in these studies limits possible conclusions 
about the appearance of strongly localized ILMs or DBs. 

A nonperturbative and qualitatively new property of 
the evolution of spatially localized excitations takes place 
if the nonlinear terms become essential. Instead of ob- 
serving a decay in space, robust and seemingly exact 
spatially localized vibrational excitations have been ob- 
served. Early results on ILMs or DBs in FPU chains have 
been obtained by Sievers, Takeno and Kisoda These 
studies as well as their followups 20] provide three types 
of observations and results, among others. 

First, they report on numerical evidence of the exis- 
tence of long lived localized excitations, with lifetimes 
much larger than the typical oscillation times 27r/cL;^. 
The internal frequencies of such ILMs are outside the 
spectrum ojq. In Figs l4lt)lSl some of these excitations are 
shown (and will be discussed in more detail below). In 



fact, as will be discussed in the next sections, these exci- 
tations are exact solutions of the equations of motion. 

Second, approximate numerical solutions were ob- 
tained using the Rotating Wave Approximation which im- 
plies that a time-periodic solution is constructed taking 
into account only its first harmonics frequency contribu- 
tion (and optionally also a dc component) and neglect- 
ing higher harmonics. The recipe is to use the ansatz 
xi{t) = Q + aicos{Q,bt)^ to insert it into the equations 
of motion (0) and to neglect all terms with higher har- 
monics which appear due to the nonlinear terms. The 
resulting set of coupled algebraic equations for the coeffi- 
cients c; , ai can be solved numerically for finite numbers 
of sites, and yields solutions similar to those observed 
in numerical simulations. This approximation is strictly 
speaking valid only for small amplitudes of the displace- 
ments. In reality however it may serve as a good estimate 
to exact solutions for rather large amplitudes as well. 

And third, for quite long simulation times moving 
ILMs have been observed, which in addition of being 
characterized by internal oscillations, propagate along 
the chain. Typically the observed propagation veloc- 
ities are smaller than the maximum group velocity 
max{\vq\) — -\/02. Also the motion of these ILMs along 
the lattice leaves excited lattice parts behind, implying 
that in the course of time these moving structures will 
slow down or disappear due to radiation of energy. 

While initially ILM excitations seemed to be connected 
to some specific properties of FPU chains, the obser- 
vation of similar localized structures in Klein-Gordon 
chains as well as in higher dimensional lattices suggested 
that the existence of ILMs is a rather generic feature for 
various anharmonic Hamiltonian lattices 0. This view 
became a well established fact due to numerous stud- 
ies during the past decade. In the following we will dis- 
cuss the existence and properties of ILMs or DBs in FPU 
chains using methods of the modern theory of localized 
excitations in discrete systems. 

DISCRETE BREATHERS IN FPU CHAINS - 
SOME DEFINITIONS 

Discrete breathers (intrinsic localized modes) are time- 
periodic spatially localized solutions of the equations of 
motion of a Hamiltonian lattice p|. In mathematical 
terms we are searching for solutions of JSJ satisfying 

xi{t + Tb)=xi{t) , pi{t + Ti,)^pi{t) , (8) 
xi^±oo ~* d± , pi^±oo ^ . (9) 

The difference d+ — d- characterizes the DB induced lat- 
tice deformation and is related to the abovementioned 
effects of thermal expansion (or contraction). Without 
loss of generality we may choose = in the following. 
The period T;, is related to the DB frequency Q.b = iii/Ti,. 
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The DB solution can be thus represented as a Fourier se- 
ries expansion with respect to time: 

xi{t)^ A.ie''^''^* . (10) 

fc— — oo 

The locahzation property @ implies 

Ak^o,i-^±oo ~* , Ak=o,i~*±oo — ^ d± . (11) 

Inserting 1)10(1 into the equations of motion (jSJ, and as- 
suming a large distance from the DB core, the lineariza- 
tion of the algebraic equations for the coefficients Aki to- 
gether with the condition (|ll|l leads to the nonresonance 
condition 

kilb ^ UJg (12) 

for all integer k. Excluding fc = for a moment, that 
condition implies > lVt^- The possible resonance for 
k — which induces nonzero static lattice deformations, 
will be discussed in more detail below. 

The stability of a DB (as for any periodic orbit) can be 
accounted for by linearizing the phase space flow around 
a given DB solution xi{t), i.e. by adding a small pertur- 
bation to it xi(t) = xi(t) + ei{t), inserting this expression 
into the equations of motion (j^J and keeping only terms 
linear in ef. 

Q = TT; , 7T; = -W"{xi - Xi-i){ei - e;_i) + 

W"{xi+i-xi){ei+i-ei) . (13) 
Equations l|13|l define a map 

(f:,')-(?oO 

which maps the phase space of perturbations onto itself 
by integrating each point over the DB period Tf,. Here 
we used the abbreviation a; = {xi, X2, xi, ...). The map 
((T^ is characterized by a symplectic Floquet matrix J^, 
whose complex eigenvalues A and eigenvectors y provide 
information about the stability of the DB. For details 
we refer to jl^. Here we note that if all eigenvalues A 
are of length one, then the DB is linearly (marginally) 
stable. Otherwise perturbations exist which will grow in 
time (typically exponentially) and correspond to a lin- 
early unstable DB. Upon changing a control parameter 
(e.g. the DB frequency) stable DBs can become unstable 
(and vice versa). Such a change of stability is appearing 
because two (or more) Floquet eigenvalues collide on the 
unit circle and depart from it. If one of the two associ- 
ated eigenvectors (or both) is spatially localized, such a 
collision and the corresponding instability are indepen- 
dent of the size of the lattice. If both eigenvectors are 
spatially delocalized, the strength of the instability de- 
pends on the size of the system and vanishes in the limit 
of an infinite system j23| . 



The absence of a cubic term (^3 (or more generally 
the case of a symmetric potential W) implies a parity 
symmetry of the interaction potential W. Consequently 
DB solutions will contain only odd harmonics in a Fourier 
expansion with respect to time, A2k,i = 0. That implies 
(i+ = (we remind the reader that d- — a,s well). For 
such a case the Floquet problem l|13() is periodic with 
period Tb/2, because the time-periodic coefficients in (fl!^ 
contain only even harmonics 2fcr2t, in a Fourier expansion 
with respect to time. 

Two eigenvectors of the Floquet matrix are corre- 
sponding to homogeneous shifts in all coordinates or like- 
wise velocities (due to the conservation of P). Their 
eigenvalues are always located at -1-1 in the complex 
plane. Two more eigenvectors of the Floquet matrix are 
corresponding to perturbations along the DB periodic 
orbit (phase mode) or along the family of DB solutions. 
Their eigenvalues are located either at +1 for asymmet- 
ric potentials W (since the Floquet mapping is performed 
for one period Tf,) or at -1 for symmetric potentials W 
(since the Floquet mapping is performed for Tb/2). 

Finally we will also use the linearized phase space flow 
dynamics around a DB solution ((13|l in order to compute 
the transmission coefficient for a small amplitude plane 
wave launched into the DB. It is a special Bloch-type 
solution of the equations H13() of the form [2^ 

00 

ei{t)= Qfce'^"'+'''''^* . (15) 

fc— — 00 

The DB solution acts as a time-periodic scattering poten- 
tial for an incoming wave with frequency ujq and gener- 
ates new channels at frequencies ujq + kftb- If + kflb 7^ 
Uqi for nonzero k then all channels are closed except for 
the open channel k — 0. Such a situation corresponds 
to elastic scattering, i.e. the energies of the incoming 
and the outgoing (transmitted and reflected) waves are 
equal. Otherwise we are confronted with inelastic scat- 
tering. The scattering setup is schematically represented 
in Figl^l For computational details the reader is referred 
to |25|. 



EXISTENCE PROOFS 

The pioneering DB existence proof of MacKay and 
Aubry [2y| was designed to apply for networks of weakly 
coupled oscillators using the implicit function theorem. 
The extension of this technique to FPU systems turned 
out to be complicated because one needs a limiting case 
where DB solutions are compact. Nevertheless it was ap- 
plied successfully to the case of a diatomic FPU chain 
with alternating heavy and light masses 27]. DB solu- 
tions were shown to exist close to the limit where the 
ratio of light to heavy masses vanishes. 
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FIG. 3: Schematic representation of the one-channel scatter- 
ing of a wave by a discrete breather. 

A special case of homogeneous potentials W with 0„i = 

t,2n where n is a given positive integer, was treated in 
to separate time and space dependence. This feature 
was then used in to provide a constructive proof of 
existence of DBs. As shown recently ,30|, this proof can 
be extended to systems with onsite potentials by adding 
Y^ii^ii^o^f + I'xf"-) to the Hamiltonian (QJ. 

An implicit proof of existence of DBs in FPU mod- 
els was provided by Aubry et al in jsj] using a varia- 
tional method. Discrete breathers are obtained as loops 
in phase space which maximize a certain average energy 
function for a fixed pseudoaction. More recently James 
[3^ proved existence of DBs in FPU models for low am- 
plitudes (energies). 

All these results provide the certainty of existence of 
DB solutions in infinite FPU models. However, due to 
the implicit character of most of these proofs a detailed 
study of DB properties has to be obtained e.g. using 
advanced computational methods. 

DB SOLUTIONS AND THEIR STABILITY 

Computational tools for studying DB properties are 
confined to the case of a finite lattice size. The typ- 
ically exponential spatial degree of localization of DBs 
yields reliable results which apply for infinite lattice size 
as well. In terms of the dynamics in phase space {xi,pi} 
we are searching for a starting point of a trajectory such 
that after integrating over a given period of time the 
trajectory reaches the starting point again. Of course 
all points on the corresponding one-dimensional man- 
ifold generated by the trajectory fulfill this condition. 
Discrete breathers are so-called isolated periodic orbits. 
That implies that fixing the total energy and mechanical 
momentum for a given DB solution, no other DB solu- 



tions are found in an infinitesimal neighbourhood of the 
DB loop. In other words, DB solutions do not belong 
to higher dimensional resonant tori in the model phase 
space. However if changes in the energy are allowed, then 
generically infinitesimal deformations of the original DB 
loop will result in a new DB loop with slightly changed 
parameters such as energy, frequency or amplitude. Slid- 
ing along such a family of DB solutions may lead to a 
stability change of the solutions. Here we use a New- 
ton method to find a DB solution at a chosen frequency 
together with open boundary conditions. For details on 
this method as well as on numerical evaluations of the 
Floquet matrix we refer to [2^ . 

The Page mode for <^3 = 

In Fig0| we show the profile of a stable DB solution 
xi{t = 0) for 02 = ^4 = 1, (l>3 = 0, = 4.5 and 80 
sites, the index I running from I = —39 to Z = 40 with 
open boundary conditions. Note that pi(t = 0) = 0. 
Because fih > o;^, the amphtudes stagger l2l|. The solu- 
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FIG. 4: Displacements xi{t = 0) for (/)2 = <?!)4 = 1, 03 = and 
Qb = 4.5 versus lattice site number I for an antisymmetric DB. 
Inset: Staggered deformation ui{t — 0) for the same solution 
on a logarithmic scale versus lattice site number I. This DB 
is stable. 

tion is antisymmetric in space, and likewise coined cen- 
tered between sites and also Page mode. In the inset we 
plot the profile of the corresponding staggered deforma- 
tion ui = {—iy{xi — xi^i) on a logarithmic scale. The 
DB is strongly localized, and only a part of the chain is 
shown. The inset shows that the DB is localized expo- 
nentially. The exponent de pen ds on the model parame- 
ters and the DB frequency |2ll. In the left part of FigO 
we show the location of the Floquet eigenvalues A in the 
complex plane. All eigenvalues reside on the unit circle. 
The two pairs of degenerate eigenvalues located at -1-1 
and -1 on the real axis (see discussion above) are shown 



6 



E 













□ 










..^011^ 


I- 






I 





-1 -0.5 0.5 1 -30 -20 -10 10 20 30 40 



Re(?i) I 

FIG. 5: Left panel: location of Floquet eigenvalues A in 
the complex plane for the DB in FigU] (crosses, diamonds, 
squares, circles). The unit circle is shown to guide the eye. 
Right panels: real part of the displacement components of the 
Floquet eigenvectors marked with the corresponding symbols 
(square and circle). 

with diamonds. Besides the Floquet continuum of ex- 
tended eigenstates (crosses) two pairs of stable spatially 
localized eigenstates are observed (shown by squares and 
circles) . In the two right panels of Fig|31 the real part 
of the displacement components of these localized eigen- 
states are plotted. 

The Sievers-Takeno mode for (/>3 = 

For the same parameters we show in FiglSjan unstable 
DB solution which is symmetric in space, or centered on 
a site and also known as the Sievers-Takeno mode. 



the Page mode, one of the pairs of localized eigenstates 
(circles) is located on the real axis off the unit circle. In 



-0.5 



-o < 


^,..---E 

> o 














"""■--c: 







-1.5 -1 -0.5 0.5 1 -30 -20 -10 10 20 30 40 



Re{X) I 

FIG. 7: Left panel: location of Floquet eigenvalues A in 
the complex plane for the DB in Fig|S] (crosses, diamonds, 
squares, circles). The unit circle is shown to guide the eye. 
Right panels: real part of the displacement components of the 
Floquet eigenvectors marked with the corresponding symbols 
(square and circle). 

the two right panels of Fig|Sl the real part of the dis- 
placement components of these localized eigenstates are 
plotted. The instability is caused by a spatially localized 
Floquet eigenvector, which deforms the Sievers-Takeno 
mode in Fig|H|in the direction of the Page mode in Fig^l 
Indeed the Sievers-Takeno mode, when perturbed in the 
full nonlinear equations of motion along the unstable Flo- 
quet eigenvector, starts to perform additional oscillations 
around the stable Page mode for times large compared 
to the DB period. 
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FIG. 6: Displacements xi{t = 0) for (l}2 — (j}4 — 1, (t>3 = and 
Qb ~ 4.5 versus lattice site number / for a symmetric DB. 
Inset: Staggered deformation ui{t = 0) for the same solution 
on a logarithmic scale versus lattice site number I. This DB 
is unstable. 

In the left part of Fig[7| we show the location of the 
Floquet eigenvalues A in the complex plane. Not all 
eigenvalues reside on the unit circle. At variance with 



The Page mode for <j)3 ^ 

li (j>3 0, then the interaction potential W becomes 
asymmetric. DB solutions will now contain in general 
all Fourier harmonics kUb including fc = 0. Depending 
on the sign of <^3 this will cause either a contraction or 
expansion of the chain. In Fig|21we show a DB solution 
which can be continued from the Page mode in Fig^up 
to 03 = 1. For that particular case d+ « —1.005. The 
DB solution can be thus viewed as a localized vibration 
which induces a nonzero kink-shaped lattice distortion. 
This DB solution turns out to be unstable. The reason 
for that is the interaction of localized Floquet eigenstates 
with the Floquet continuum (the eigenstates which cor- 
respond to plane waves far away from the DB). Such 
oscillatory instabilities can be characterized by the wave 
numbers Qos which correspond to the extended Floquet 
state interacting with the originally localized one [2^ . In 
the left part of Fig|51we show the location of the Floquet 
eigenvalues A in the complex plane. Not all eigenval- 
ues reside on the unit circle. The two pairs of degener- 
ate eigenvalues, which were located at -f 1 and -1 on the 
real axis (see discussion above) for 03 = 0, are now all 
located at -1-1 and shown with diamonds. Besides the 
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FIG. 8: Displacements xi{t — 0) for 02 = 03 = 04 = 1 
and Qb = 4.5 versus lattice site number /. Inset; Staggered 
deformation ui{t = 0) for the same solution on a logarithmic 
scale versus lattice site number I. This DB is unstable. 



position and starts to move along the lattice. 

The Sievers-Takeno mode for 03 7^ 

When continuing the Sievers-Takeno mode from 03 = 
to nonzero values of 03, the strong instability caused 
by a localized Floquet eigenstate for 03 = remains. It 
also acquires the above discussed oscillatory instabilities 
of the Page mode. At the same time the spatial profile 
of the mode ceases to show up with any symmetry. 

DB stability properties for different frequencies 

In the following we present results on the stability of 
the Page mode DB upon variation of the DB frequency 
rif). The Sievers-Takeno mode DB yields similar results, 
except for keeping its strong localized instability dis- 
cussed above. 
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FIG. 9: Left panel: location of Floquet eigenvalues A in 
the complex plane for the DB in Fig|H| (crosses, diamonds, 
squares, circles). The unit circle is shown to guide the eye. 
Right panels: real part of the displacement components of the 
Floquet eigenvectors marked with the corresponding symbols 
(square and circle). 

Floquet continuum of extended eigenstates two pairs of 
eigenstates (which were localized for the Page mode with 
03 — 0) are observed to strongly interact with the Flo- 
quet continuum (shown by squares and circles). While 
the small deviation of the corresponding eigenvalues from 
the unit circle is of the order of 0.01. ..0.03, we checked 
that it is independent of the system size by increasing 
the number of sites from 80 to 640 (see also [23j). In the 
two right panels of FigElthe real part of the displacement 
components of these resonating and unstable eigenstates 
are plotted. 

When the Page mode is perturbed along the oscillatory 
instability eigenvector in the full nonlinear equations of 
motion, it deforms initially very slowly compared to the 
DB period, due to the weak amplitude of the instability. 
After reaching a critical threshold in the perturbation 
amplitude, the mode depins from its original lattice site 
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FIG. 10: All Floquet phases 6 and the squared eigenvalue 
length |Ap of the unstable Floquet eigenvectors for a DB with 
02 = 04 = 1 and 03 = 0.5 as a function of the DB frequency 
Qb. Eigenvalues which correspond to finite size instabilities 
are shown using grey symbols. 

In Fig^l we plot the dependence of the arguments 9 
and the squared absolute values |Ap of the irreducible 
part of the Floquet matrix eigenvalue spectrum A = 
|A|exp(z6') as a function of the DB frequency fif, for 
02 = 04 = 1 and 03 = 0.5. These results are ob- 
tained for a DB family which is continued from the Page 
mode in Fig^ We first note that lUt^ = 2 and for 
LUj^ < flh < 2uJt^ the arguments 9 cover the whole unit 
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circle interval. The corresponding interactions between 
different extended Floquet eigenstates lead to finite size 
instabilities [l^jjlj which are marked with grey sym- 
bols in the lower panel of Fig^| These instabilities will 
disappear for an infinite system, i.e. the corresponding 
values |A| — 1 ~ l/N where N is the number of sites [isj . 
For that DB frequency interval scattering of waves will 
be inelastic for wavenumbers q, q' satisfying ojq + Q.b — ujq> 
|23 |. Note that this inelastic two-channel scattering will 
survive in the limit of an infinite system, opposite to the 
finite size instabilities of the DB itself. 

The mentioned oscillatory instabilities induced by the 
interaction between two pairs of originally localized 
eigenstates and the continuum are clearly observable, and 
cause deviations of the corresponding absolute values |A| 
from unity. As conjectured in Jsll, these instabilities sur- 
vive in the limit of an infinite lattice. The first pair of 
these eigenstates exits the continuum at fib w 4.6. The 
corresponding eigenvalues return to the unit circle and 
approach —1 in the complex plane at Qb ~ 12, collide 
with each other, separate on the real axis causing an- 
other instability (cf. the corresponding ellipse structure 
in Fig llU|) and merge again at -1 for fif, w 14. The second 
pair exits the continuum at fli, w 11.5. For sufficiently 
large values of > 14 the DB becomes stable. 

The discussed instabilities for the Page mode are 
solely caused by the asymmetric term 03 . Indeed, for 
03 = the abovementioned periodicity of H13(l with 
Tfc/2 excludes the possibility of finite size instabilities or 
equally inelastic multichannel scattering, since the condi- 
tion ujq + 2flb — ujq' can not be satisfied due to ^b > ij->-n- 
The oscillatory instabilities caused by the interaction of 




2 4 6 8 10 12 14 16 18 20 

FIG. 11: All Floquet phases 9 for a DB with 02 = 04 = 1 
and 03 = as a function of the DB frequency Q.i,. Note that 
the Floquet map is obtained here by integrating 113^ over the 
time Tb/2. 

localized Floquet eigenstates with the Floquet continuum 
are removed for the same reason. The dependence of the 



eigenvalue phases Q on Vlb for a Floquet matrix defined 
by integrating (|13|l over Tb for 03 = is practically iden- 
tical to the case 03 = 0.5 (upper panel in Fig lTUI) . How- 
ever the correct irreducible Floquet matrix for 03 = is 
obtained by integrating ((T^ over the true period Tb/2. 
The phase variation shown in Fig^J demonstrates that 
the localized Floquet states are indeed separated from 
the Floquet continuum. 

RESONANT WAVE SCATTERING BY DBS 

As mentioned above, the linearized phase space flow 
equations (|13|l give also information about the complex 
transmission amplitude tq and the transmission coeffi- 
cient T{q) = While the numerical technique and 
also a number of results concerning wave scatterin g by 
FPU breathers have been reported recently |2^.[25|.|33|. 
we want to address here one aspect which was not yet 
discussed. 

We show in Fig^jthe g-dependence of T{q) for three 
different frequencies fif, = 4.5,4.7, 10 of the Page mode 
DB with 02 = 04 = 1 and 03 = 0.5, for which we gave the 
Floquet spectral data in Fig^] The perfect transmission 
T{q ~ {)) = 1 for all cases is due to the conservation of 
total mechanical momentum |25). Equally the vanishing 
of the transmission T{q = tt) = is due to the vanishing 
of the group velocity {T)) at the band edge q = ir [23| . 
However the resonant perfect transmission and reflection 
peaks observed for g ^ 0, tt are due to Fano resonances 
|25l|.[33|. We remind the reader that Fano resonances are 
induced by localized states interacting with a continuum 
of scattering states [^. 

This is precisely the situation observed in the case of 
an asymmetric FPU potential. For 03 = the Page 
mode DB has a localized Floquet eigenstate well sepa- 
rated from the Floquet continuum. However any nonzero 
value of 03 induces an interaction between this localized 
state and the continuum mediated by the DB. This in- 
teraction causes an oscillatory instability in the Floquet 
matrix spectrum. The dashed vertical lines in all three 
panels in Fig^|indicate the location of the oscillatory in- 
stability qos as observed in Fig^| In all three cases the 
value of qos is very close to the locations of perfect trans- 
mission and reflection. In Fig^Jwe plot the g- values of 
the oscillatory instability of the DB qos from Fig^] to- 
gether with the positions of the perfect transmission and 
reflection as functions of fib. The obtained correlation 
between these resonant g- values is evident. 

From a general point of view, the DB acts as a time- 
periodic scattering potential for incoming waves. The 
time-periodic scattering potential generates locally new 
frequencies from an incoming wave which correspond to 
closed channels. These closed channels provide addi- 
tional propagation ways for the incoming wave. In that 
sense a time-periodic scattering potential in a strictly 
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FIG. 12: Transmission coefficient T versus wave number q of 
an incident wave for (j>2 = 4'i = ^ and 03 = 0.5 scattered by a 
Page mode DB. From top to bottom: Q.b = 4.5, 4.7, 10. Tfie 
vertical dashed line indicates the position of the oscillatory 
instabilities qoe- 



one-dimensional system acts similar to a local enlarge- 
ment of the system dimensionality around the scatterer 
[i^l,!!^. Such additional propagation channels allow for 
wave interference, and may also cause a total destructive 
interference leading to a total reflection of the incoming 
wave, which is another way of understanding the origin 
of the observed Fano resonance [s^l . Note that in a strict 
sense the correlation between an oscillatory instability of 




FIG. 13: The dependence of the position of an oscillatory in- 
stability Qos (diamonds) and the associated locations of per- 
fect transmission (solid line) and perfect reflection (dashed 
line) on Q,b for a Page mode DB with 02 = 04 = 1 and 
03 = 0.5. 



a DB solution and a corresponding Fano resonance holds 
only for the limit of weak coupling between the local- 
ized Floquet state and the Floquet continuum . And 
only in that limiting case we can expect a nearby loca- 
tion of perfect reflection (Fano resonance) and perfect 
transmission. The case of strong interaction can lead to 
the disappearance of an oscillatory instability in the Flo- 
quet spectrum and of the perfect transmission peak, but 
simultaneously to a remaining of the Fano resonance of 
perfect reflection [s^ ]. 



THE FREQUENCY LIMITS OF THE DB 
SOLUTION FAMILIES 



When the frequency fib is varied along a DB family, we 
may consider the two limiting cases fifc ^ oo and — > 
dj^r. The high frequency limit implies large amplitudes 
inside the DB core. The leading order contribution in 
the equations of motion will then be due to the 04 term. 
The core profile of these high frequency DBs will be very 
similar to the DB cores in Figs l4l6l reefardless the values of 
02 and 03 which will only affect the tail characteristics 
[13; [1^7 [13' Opposite to that, the low frequency 
limit is depending on the value of 03 [ssf.fssj. 

In Fig^^we plot the dependence of the DB energy ^ 
on its frequency for 03 = 1 and 03 = 0.5 and 80 sites. 
For 03 = 0.5 the energy tends to zero as Jlfc tends to 
cjjr. However for 03 = 1 the DB energy passes through a 
minimum as the frequency is lowered, and increases a gain 
with further approaching the limiting value a;,r . |38l| . 
These results are intimately connected with the modu- 
lational instability of plane waves at small amplitudes 
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FIG. 14: DB energy E versus DB frequency Qi, for 02 ~ 
(l>4 — 1 and <j!>3 = 0.5 (lower curve) and (j>3 = 1 (upper curve). 



Insets: Displacements xi{t = 0) for ilj, 
site number I. Upper left inset for (ps 
for <^3 = 0.5. 



2.001 versus lattice 
1, lower right inset 



[13], 15], 16"!. Assuming that a DB family has a zero am- 
phtude hmit, it follows that in a one-dimensional chain 
as discussed here its energy will also vanish in that limit 
[395. It has been conjectured (and later confirmed 
numerically [isj ') that a small-amplitude DB occurs due 
to an instability (bifurcation) of the q = n plane wave. 
The rigorous analysis of this plane wave bifurcation gives 
an inequahty (13t|.[l5 | .[l ^ .|3 ^ 



(16) 



which has to hold in order to obtain the modulational 
instability of the q = ir mode and by conjecture a low- 
amplitude DB solution. Since inequality (|16|l is satisfied 
for (jj^ = (f'i — ^ (jjT, = 0.5, we do observe the low 
amplitude (and thus low energy) limit of the DB family 
in Figll4l However choosing 412 — fpi — ^ and 03 = 1, 
the inequality (|16|l is violated. Consequently the q — tt 
mode does not experience a low amplitude modulational 
instability, and low amplitude (and low energy) DBs do 
not exist. Nevertheless, as explained above, the high fre- 
quency limit of such an FPU model will always allow for 
DB solutions of a well known asymptotic shape. Conse- 
quently the energy versus frequency curve in Fig 1 141 has 
to show up with a minimum. It is instructive to compute 
the DB profiles for 11;, close to lUt^ for the two cases. These 
profiles are shown in the insets in Fig^l For ^3 = 0.5 
the low energy DB delocalizes and takes the shape of a 
g = TT plane wave (or better to say a standing wave, due 
to the open boundary conditions). The lattice deforma- 
tion d+ « —0.03 tends to zero when increasing the system 
size and further approaching 2. For 03 = 1 the DB 

with frequency close to has still large amplitudes, in- 
cluding a nonzero and well pronounced lattice distortion 



d+ ~ —1. While both DB branches have weak oscilla- 
tory instabilities as shown in Fig llOl the dashed part of 
the upper curve in FigEl indicates that DBs from that 
part of the branch experience an additional instability 
with corresponding eigenvalues located on the real axis. 
When decreasing 03 and approaching 4>'3^^ — \/30204/2, 
which changes the sign of the l.h.s. of inequality Hlt)|l . 
the value of d+ for DBs with frequencies close to the 



linear spectrum (fib - 
threshold disappears. 



2) tends to zero and the energy 



DISCUSSION 

The results of this work can be continued - and have 
been continued - in many different directions. The unsta- 
ble Sievers-Takeno mode DB, when properly perturbed, 
lead to DB-like excitations which propagate along the 
lattice Some perturbative analytical results give 

a partial understanding of this effect [4l| . Still there is 
strong evidence of the nonexistence of exact and spatially 
localized moving DBs |4^ . 

The inequality H16|l can be used to predict the existence 
or noexistence of DB solutions for various other FPU-like 
models. A consequence is that DB excitations are not 
expected to exist for the Toda chain. A similar conclusion 
can be obtained for a so-called Roto- FPU model with 
W{x) = 1 — cosx. However such a system of coupled 
rotors allows for rotobreathers, i.e. solutions where a few 
degrees of freedom are in a rotating state, while the rest 
of the system is in a spatially localized oscillating state 
M (cf. FigCttb)). 

Discrete breathers are not an exclusive property of one- 
dimensional lattices. They exist equally in higher dimen- 
sional lattices Q , where they show up with nonzero lower 
energy thresholds js^ , independent on further model pa- 
rameters. This is true also for generalizations of FPU 
models to higher lattice dimensions. An interesting result 
concerns the case of asymmetric interaction potentials W 
and higher-dimensional FPU lattices. As shown in [i^ ]. 
the requirement to have a finite energy for a DB leads to 
a static lattice deformation with a dipole symmetry and 



a spatial decay law 



where r is the distance from 



the breather core and d the lattice dimension. 

Finally we mention a number of studies of FPU 
breathers emerging from perturbed extended states via 
modulational instability or induced by fluctuations in 
thermal equilibrium [l3j|4a|- These studies show clearly 
that DB solutions are not only interesting mathematical 
objects, but essential in order to understand and describe 
the properties of nonlinear lattice dynamics in thermal 
equilibrium and during relaxational processes. 
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